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I .  INTRODUCTION 


A  loss  of  coolant  accident  (10CA)  is  defined  as  a  failure  of  the 
reactor  coolant  system  resulting  in  a  partial  or  complete  release  of 
coolant  to  the  containment  vessel,  thereby  interrupting  the  normal  mech¬ 
anism  for  removing  heat  from  the  reactor  core.  If  emergency  coolant  were 
not  provided,  clad  and  fuel  melting  would  eventually  occur  as  a  result 
of  residual  heat  generation  in  the  fuel  and  possible  metal  water  and 
other  chemical  reactions.  The  purpose  of  the  Emergency  Core  Cooling 
System  (ECCS)  is  to  provide  adequate  core  cooling  for  the  period  imme¬ 
diately  following  a  break  and  to  assure  that  the  core  is  recovered 
and  remains  covered  so  that  the  residual  heat  can  be  removed. 

The  present  MH-  1A  Emergency  Core  Cooling  System  consists  of  an  emergency 
coolant  injection  tank  and  two  high  pressure  coolant  charging  pumps. 

Emergency  operation  of  the  injection  tank  and  one  coolant  charging  pump  is 
initiated  by  a  coincident  signal  of  low  pressurizer  pressure  and  high 
containment  building  pressure.  Cooling  water  from  the  injection  tank  is 
delivered  through  a  check  valve  to  the  reactor  vessel  downcomer  region 
when  the  pressure  in  the  inlet  plenum  drops  below  approximately  240  psia. 
Driving  head  for  this  injection  consists  of  a  regulated  gas  pressure 
within  the  injection  tank  plus  the  elevation  head.  The  injection  tank  is 
designed  to  recover  the  core  by  bottom  flooding  within  approximately  ten 
minutes.  The  coolant  charging  pumps  serve  to  make  up  coolant  lost  due  to 
a  small  leak  or  as  a  result  of  boil-off  after  recovering  of  the  core  and 
provide  long  term  emergency  cooling  to  remove  decay  heat.  Table  1  summarizes 
the  characteristics  of  these  systems  which  are  described  in  more  detail 
in  Appendix  A  and  References  3  and  5. 

This  report  evaluates  the  effectiveness  of  the  present  MH-1A  ECCS 
during  a  postulated  loss  of  coolant  accident  and  presents  information  to 
establish  the  required  performance  of  a  system  meeting  the  AEC's  "Interim 
Acceptance  Criteria  for  Emergency  Core  Cooling  Systems"  (Reference  1)  published 
in  the  Federal  Register  June  29,  1971.  All  calculations  were  performed  in 
accordance  with  the  AEC  evaluation  model  Appendix  A,  Part  1  of  the  above 
reference  using  the  recommended  compute r  codes  RELAP-3  (Reference  7)  for 
reactor  blowdown  calculations  and  THETA  1-B  (Reference  8)  for  the  fuel 
element  heat  up  transients.  Analyses  reported  herein  cover  break  sizes 
from  0.18  ft2  up  to  the  double  ended  severance  of  the  largest  reactor 
coolant  pipe.  The  analyses  were  performed  at  110  percent  of  rated  power 
(49.5  MWt)  and  at  a  peak  linear  heat  rate  of  21.1  kw/ft. 
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TABLE  1 


SUMMARY  OF  MH-1A  EMERGENCY  CORE  COOLING  SYSTEM  CHARACTERISTICS 


Coolant  Charging  Pumps 


Number  of  pumps' 

Capacity  per  pun?) 

Actuation  (coincident) 

Emergency  power  source 

Startup  delay  for  diesel 
Point  of  Injection 
Coolant  source 
Capacity  of  tank 


Emergency  Injection  System 

Number  of  injection  lines 
Emergency  injection  tanks 
Method  of  injection 
Operating  pressure 
Injection  tank  capacity 
Minimum  emergency  coolant 
Coolant  temperature 
Check  valves  in  injection 
line 

Isolation  valves  in  injection 
line 

Isolation  valves  in  nitrogen 
pressurizing  line 
Method  of  actuation 
(coincident) 

Design  injection  rate 
Point  of  injection 


2 

8.5  gpm 

High  coat,  pressure  (1/2) 

Low  pressurizer  press  (1/2) 
Diesel  generator 
^  30  sec 

Reactor  Vessel  downcomer 
Emergency  water  storage  tank 
10,000  gallons 


1 

1 

pressurized  gas 
240  psia 
2600  gal 
1000  gal 
85°F 

L 

2 

2 

{High  cont.  pressure  (1/2) 
Low  Pressurizer  press  (1/2) 
^  100  gpm 

Reactor  Vessel  downcomer 
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II-  SUMMARY  AND  CONCLUSIONS 


The  performance  of  the  ECCS  presently  installed  at  the  MH-1A  has 
been  analyzed  for  large  and  intermediate  size  breaks  within  the  reactor 
primary  system  using  the  methods  and  assumptions  specified  in  the  AEG 
Interim  Policy  Statement  on  Emergency  Core  Cooling  Systems  (Reference  1). 

The  parameters  used  in  the  analyses  are  consistent  with  the  data 
reported  in  the  MH-1A  SAR  (Reference  3)  and  MH-2A  Core  4  Nuclear  Analysis 
(Reference  9).  Table  2  summarizes  the  MH-1A  thermal  parameters  related 
to  the  LOCA.  The  computer  programs  and  methods  used  in  analyzing  the 
effectiveness  of  the  ECCS  are  discussed  in  Section  IV.  The  major  input 
parameters  and  additional  data  developed  for  the  LOCA  analysis  are  reported 
in  the  Appendices  A  through  H  and  the  results  of  the  analyses  for  the 
present  system  are  given  in  Section  V.  Based  on  these  results  it  is 
concluded  that  the  criteria  identified  in  Section  IV-A  of  the  AEC  Interim 
Policy  Statement  are  not  met  by  the  present  MH-1A  ECCS.  Further  parametric 
analyses  were  performed  to  establish  the  basic  performance  criteria  neces¬ 
sary  for  the  Improvement  of  the  ECCS  to  meet  the  AEC  interim  criteria. 

The  main  conclusions  regarding  the  improvement  of  the  MH-1A  ECCS  as  a 
result  of  these  analyses  are  discussed  in  Section  VI. 

TABLE  2 

SlflMAKY  OF  MH-1A  THERMAL  DESIGN  PARAMETERS 


Thermal  power  rating,  MWt 
Number  of  fuel  assemblies 
Fuel  rods  per  assembly 
Cladding  outside  diameter,  in 
Fuel  pellet  diameter,  in 
Fuel  rod  pitch,  in 
Cladlng  material 
Average  fuel  specific  power 
kw/kg  of  UO2 
kw/ft  of  rod 

Adiabatic  heating  rate  of  fuel 

at  average  specific  power  °F/sec 
Average  rod  power,  kw 
Peak  to  average  beat  flux  ratios 
Axial 
Radial 
Total 

Average  fuel  temperature  at 

highest  heat  flux  position,  °F 
Core  average  fuel  temperature  at 
rated  power,  °F 


45 

32 

104 

.507 

.4565 

.654 

348  stainless  steel 

13.6 

4.35 


66.5 

13.1 

1.98 

1.95 

3.86 

2510 

1060 
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III.  CONFORMANCE  WITH  ARC  INTERIM  Ml. ICY  SiXilSttl 


A.  GENERAL 


The  analysis  of  the  KH-lA  LOCA  was  perforaed  in  accordance  with  the 
AEC  Policy  Statement  for  eaergency  core  cooling  systeas  (Reference  1). 

The  criteria  for  the  acceptance  of  the  ECCS  as  contained  in  the  AEC 
policy  stateaent  are  set  if  the  results  of  the  LOCA  are  limited  as  follows: 

"1.  The  calculated  maximum  fuel  eleaent  cladding  temperature  does 
not  exceed  2300° F.  This  limit  has  been  cho6en  on  the  basis  of 
available  data  for  enbrittlement  and  possible  subsequent  shattering 
of  the  cladding.  The  results  of  further  detailed  experiments 
could  be  the  basis  for  future  revision  of  this  limit. 

2.  The  amount  of  fuel  eleaent  cladding  that  reacts  chemically 
with  water  or  steal  does  not  exceed  1  percent  of  the  total  amount 
of  cladding  in  the  reactor. 

3.  The  clad  temperature  transient  is  terminated  at  a  time  when 
the  core  geometry  is  still  amenable  to  cooling,  and  before  the 
cladding  is  so  enfcrittled  as  to  fail  during  or  after  quenching. 

4.  The  core  temperature  is  reduced  and  decay  heat  is  removed  for 
an  extended  period  of  time,  as  required  by  the  long-lived  radio¬ 
activity  remaining  in  the  core." 

This  section  discusses  the  application  of  the  AEC  evaluation  model 
contained  in  Appendix  A,  Part  I  of  the  interim  policy  statement  to  the 
analysis  of  MH-1A  ECCS  and  summarizes  the  analytical  procedures  and 
assumptions  used  in  the  evaluation  model. 

The  analyses  were  performed  for  the  entire  break  spectrum,  from  0.18 
ft2  up  to  and  including  the  double-ended  severance  of  the  largest  pipe 
of  the  reactor  coolant  pressure  boundary.  Various  break  sizes  and  locations 
were  analyzed  to  identify  the  most  severe  case.  Table  3  summarizes  the 
spectrum  of  breaks  analyzed.  Evaluation  of  the  present  MH-1A  ECCS  in 
terms  of  the  single  failure  criterion  is  discussed  in  Reference  5. 

In  accordance  with  the  AEC  evaluation  model  the  analyses  were  performed 
in  three  parts  consisting  of  system  blowdown,  refill  and  reflooding,  and 
core  heat  up  phases.  The  thermal  hydraulic  calculations  during  blowdown 
were  performed  using  the  RELAP-3  computer  program.  Four  nodes  were  used 
to  describe  the  reactor  core,  nine  nodes  for  the  primary  side  of  the  steam 
generator  and  one  node  for  the  containment.  The  re flood  transient  was 
performed  using  a  series  of  engineering  calculations  and  the  fuel  element 
heat  up  calculations  were  performed  using  the  THETA  1-B  computer  code.  Six 
radial  fuel  nodes,  two  radial  cladding  nodes  and  eleven  axial  fluid  nodes 
were  used. 
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TABLE  3 


BREAK  SPECTRUM  FOR  MH-1A  LOSS  OF  COOLANT  ANALYSIS 

1.  Double-ended  hot  leg  break  (guillotine) 

2.  Double-ended  hot  leg  break  (slot) 

3.  Single-ended  hot  leg  break  (slot) 

4.  Double-ended  break  at  the  pump  discharge  (guillotine) 

5.  Double-ended  break  at  the  pu^>  discharge  (slct) 

6.  0.6  Double-ended  break  at  the  pump  discharge  (slot) 

7.  0.6  Double-ended  break  at  the  piap  inlet  (slot) 

8.  0.5  ft2  Break  at  the  pimp  discharge  (slot) 

9.  0.18  ft2  Pressurizer  surge  line  break 
B.  ASSUMPTIONS  AND  PROCEDURES 


The  assumptions  and  procedures  specified  with  respect  to  the  AEC 
evaluation  model  axe  repeated  below  along  with  the  section  in  this  report 
In  which  their  application  to  the  MH-1A  analysis  is  discussed: 

"1.  Pump  model  -  The  pump  resistance,  K,  used  for  analysis  should  be 
fully  justified.  The  effect  of  pump  speed  upon  K 
should  be  considered.  The  more  conservative  of 
two  assumptions  (locked  or  running)  should  be  used 
for  the  pump  during  the  blowdown  calculation  (Section 
IV-B-4) . 

"2.  Break  characteristics  -  For  large  breaks  in  the  range  0.6  to  1.0 
times  the  total  area  of  the  double-ended  break  of 
the  largest  cold-leg  pipe,  two  break  models  should 
be  used.  The  first  model  should  be  the  double- 
ended  severance  (guillotine),  which  assumes  that 
there  is  break  flow  from  both  ends  of  the  broken 
pipe,  but  no  conmun'.cation  between  the  broken  ends. 
The  second  model  should  assume  discharge  from  a 
single  node  (split)  (Section  IV-B-2). 
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■'3.  A  break  discharge  coefficient  (Cq)  of  1.0  should  be  used  for 
all  break  sizes  (Section  IV-B-2). 

"4.  Decay  heat  -  The  decay  heat  curve  described  in  the  proposed  ANS 

Standard,  with  a  20  percent  allowance  for  uncertainty, 
should  be  used.  The  fraction  of  decay  heat  generated 
in  the  hot  rod  should  be  considered  to  be  100  percent 
of  this  value  unless  a  smaller  value  Is  justified 
(Section  IV-D-5). 

"5.  Time  to  departure  from  nucleate  boiling  -  use  any  calculated 
option  in  the  code  (Section  IV-D-3). 

"6.  Heat  transfer  after  departure  from  nucleate  boiling  -  use  pro¬ 
grammed  transition  boiling  correlation  option 
(Section  IV-D-3) . 

"7.  Film  boiling  heat  transfer  -  use  Groeneveld  correlation  (equation 
5.7  of  AECL  -  3281,  December  1969  (Section  IV-D-3). 

"8.  Metal-vater  reaction  rate  -  use  the  Baker-Just  equation,  with  a 
coefficient  of  1.0  (Section  IV-D-8). 

"9.  Core  flow  -  use  0.8  x  RELAP  smoothed  flow  at  the  junction  which 
is  entering  core.  If  flows  are  opposed,  use  zero 
flow  (Section  IV-D-4). 

"10.  Enthalpy  and  pressure  -  use  entering  plenum  conditions  (Section 
IV-D-4) . 

"11.  Accumulator  bypass  -  for  cold  leg  breaks,  all  of  the  water  injected 
by  the  accumulators  prior  to  end-of-blowdown  shall 
be  assumed  to  be  lost.  In  this  context  the  end-of- 
blowdown  shall  be  specified  as  the  time  at  which 
zero  break  flow  is  first  computed  (Section  IV-B-6). 

"12.  Reflood  -  A  calculation  for  the  reflooding  heat  transfer  should 

be  performed.  The  containment  back  pressure  assumed 
for  the  analysis  should  not  be  higher  than  the 
initial  pre-break  pressure  plus  80  percent  of  the 
increase  in  pressure  calculated  for  the  accident 
(Section  IV-C).  The  following  items  should  be 
constraints  on  the  calculations: 

"a.  No  steam  flow  should  be  permitted  in  intact  loops  during  the 
time  period  that  accumulators  are  Injecting. 
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Core  exit  quality  slould  be  calculated  from  entering  mass 
flow  rate  and  nominal  FLECHT  heat  transfer. 

Pump  resistance,  K,  should  be  calculated  on  the  basis  of  a 
locked  rotor. 

The  effects  of  the  nitrogen  gas  in  the  accumulator,  which  is 
discharged  following  accumulator  water  discharge,  should  be 
taken  into  account  in  calculating  steam  flow  as  a  function 
of  time. 

The  pressure  drop  in  the  steam  generator  should  be  calculated 
with  the  existing  fluid  conditions  and  associated  loss 
coefficients . 

All  effects  of  cold  injection  water,  in  either  a  hot  or  cold 
leg,  on  steam  flow  and  AP,  should  be  included  in  the  calculation. 

The  heat  transfer  coefficient  during  reflood  should  be  derived 
from  FLECHT  data." 
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IV.  METHOD  OK  ANALYSIS 


A.  GENERAL 


The  major  events  occurring  during  a  LOCA  are  (1)  system  depressurization, 
(2)  emergency  coolant  injection  and  core  reflooding  and  (3)  fuel  and  clad¬ 
ding  heat  up  until  core  cooling  becomes  effective.  The  analysis  considers 
each  of  these  three  aspects  during  the  accident.  First,  the  blowdown  cal¬ 
culation  determines  the  primary  system  hydraulic  response  and  core  power 
transient  during  system  depressurization.  This  analysis  provides  input 
for  the  calculation  of  the  core  cooling  during  blowdown  and  basic  informa¬ 
tion  regarding  the  operation  of  the  safety  injection  system.  The  blowdown 
analysis  is  performed  with  the  RELAP-3  computer  program.  A  re  flooding 
analysis  is  initiated  at  the  end  of  blowdown  which  considers  rate  of 
coolant  injection,  core  reflooding  rate,  and  heat  transfer  rates  following 
blowdown.  The  core  temperature  transient  is  then  analyzed  with  THETA-1B 
to  determine  the  cladding  temperature  response  and  the  extent  of  metal- 
water  reaction.  Input  for  this  analysis  is  determined  from  the  blowdown 
and  reflooding  analyses. 

B.  BLOWDOWN  CALCULATION 


The  RELAP-3  computer  program  (Reference  7)  is  used  to  calculate  the 
behavior  of  the  primary  system  during  subcooled  and  saturated  blowdown. 

The  code  calculates  transient  pressures,  flow  rates,  enthalpies  and  mass 
inventories  throughout  the  primary  system,  and  determines  the  core  power 
transient  and  energy  removal  rates  during  system  depressurization.  The 
code  employs  a  one  dimensional  analysis  of  the  entire  reactor  loop  re¬ 
quiring  a  geometric  description  in  terms  of  a  series  of  volumes  and  junctions. 
Pump  characteristics,  core  and  steam  generator  heat  transfer,  reactor 
kinetics,  and  emergency  injection  are  included  in  the  model. 

Before  the  break  occurs,  the  plant  is  assumed  to  be  operating  in  an 
equilibrium  condition  with  a  system  pressure  of  1500  psia,  a  reactor  coolant 
outlet  temperature  of  532°F,  .and  a  primary  coolant  flow  rate  of  3.88  x  106 
lb /hr  (9750  gpm).  A  reactor  power  level  of  110  percent  of  rated  power  is 
assumed  for  all  cases. 

1.  Primary  System  Representation.  In  the  RELAP-3  model  for  the  MH-1A, 
the  primary  system  is  divided  into  nodes  with  flow  paths  connecting  the 
nodes.  The  reactor  vessel  downcomer  region  is  represented  by  an  annulus 
extending  from  the  upper  plenum  region  to  the  bottom  of  the  inlet  plenum. 

The  lower  plenum  connects  to  the  downcomer  region  and  extends  up  to  the 
bottom  of  the  active  fuel.  The  active  core,  region  is  divided  into  four 
nodes  of  equal  volume.  Each  volume  is  initialized  at  the  proper  fluid 
temperature  based  on  the  assumed  power  distribution.  The  power  input  to 
each  of  the  axial  core  regions  is  based  on  the  relative  power  generated 
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in  chat  region  and  Che  segment  length.  The  upper  plenum  and  hot  leg  are 
each  represented  by  separate  nodes.  The  primary  side  of  the  steam  genera¬ 
tor  is  represented  by  nine  nodes;  one  node  represents  the  steam  generator 
inlet  plenum  and  another  the  outlet  plenum.  Seven  nodes  represent  the 
active  heat  transfer  portion  of  the  steam  generator  tubes.  The  cold  leg 
is  represented  by  a  node  in  which  the  primary  coolant  pumps  are  located. 
Additional  nodes  are  assigned  for  the  pressurizer,  pressurizer  surge  line, 
and  containment.  A  fill  junction  is  provided  in  the  reactor  vessel  down¬ 
comer  to  simulate  emergency  coolant  injection.  Figure  1  shows  the  repre¬ 
sentation  for  the  single  ended  hot  leg  break.  Other  models  are  similar 
except  for  break  type  and  location.  Further  details  of  the  primary  sys¬ 
tem  representation  are  discussed  in  Appendix  A. 

2.  Break  Characteristics.  The  limiting  mass  flow  through  the  break 
junction  in  RELAP-3  is  defined  by  Moody's  two-phase  choked  flow  model 
(Reference  10).  The  flew  through  the  junction  is  chosen  as  the  smaller 
of  the  inertial  flow  or  the  choked  flow.  For  choking  to  occur  the  con¬ 
tainment  pressure  must  be  less  than  the  throat  pressure.  A  discharge 
coefficient  of  1.0  is  used  in  RELAP-3  for  all  break  sizes. 

Two  break  models  were  used  for  the  analysis  of  large  breaks.  The 
first  model  was  a  double-ended  severance  of  the  pipe  (guillotine),  which 
assumed  that  there  was  break  flow  from  both  ends  of  the  pipe  but  no 
communication  between  the  broken  ends.  This  type  of  break  was  simulated 
by  assuming  leakage  from  the  adjacent  nodes  with  a  valve  in  the  junction 
joining  them.  This  valve  is  closed  at  the  time  the  break  occurs.  The 
second  model  assumed  discharge  from  a  single  node  (slot). 

3.  Reactor  Trip.  Reactor  trip  setpoints  with  appropriate  instrument 
delays  to  the  beginning  of  rod  motion  can  be  included  in  the  RELAP-3 
model.  For  the  LOCA,  reactor  trip  was  assumed  to  be  initiated  by  low 
pressurizer-  pressure  signal  (1250  psia  including  error)  with  a  .500  second 
delay  to  rod  insertion.  However,  reactor  shutdown  due  to  voids  occurred 
in  general  before  control  rod  motion  was  initiated,  for  the  break  sizes 
analyzed. 

4.  Pump  Representation.  The  pump  head  versus  flow  curve  for  the 
primary  coolant  pumps  input  into  RELAP-3  was  established  from  the  manufac¬ 
turer's  test  data  on  the  as-built  pumps.  The  test  data  was  used  to  develop 
a  nominal  speed  characteristic  curve  for  the  range  of  flow  encountered 

in  the  LOCA.  At  flow  rates  greater  than  the  flow  corresponding  to  zero 
pump  head,  the  pump  head  was  conservatively  assumed  to  be  zero.  Flow 
coastdown  following  a  pump  trip  was  approximated  with  a  pump  coastdown 
multiplier  consistent  with  the  loss  of  flow  transient  reported  in  Reference 
11.  Pump  input  data  is  discussed  in  Appendix  C. 
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The  effect  of  puny  assucytions  on  i he  core  neat  removal  rate 
during  bio /down  was  evaluated.  The  most  conservative  assumptions  are 
those  whica  result  in  the  least  pump  resistance  and  hence  the  most  rapid 
blowdown.  For  all  cases  except  the  cold  leg  break  at  the  pump  inlet, 
the  pumps  were  assumed  to  remain  running  with  the  pump  head  following 
the  characteristic  curve  until  cavitation  occurred.  Zero  puny  head 
occurred  after  cavitation.  For  the  break  at  the  pump  inlet  the  pumps 
were  assumed  to  trip  at  the  time  the  break  occurred  due  to  loss  of  off¬ 
site  power.  After  coastdown  the  resistance  corresponding  to  a  free  spinning 
rotor  was  assumed.  Conyarison  of  results  indicates  that  the  assumption  of 
pumps  running  is  conservative,  resulting  in  higher  fuel  temperatures  at 
the  end  of  the  blowdown. 

5.  Reactor  Kinetics.  Reactor  power  during  blowdown  is  determined 
from  the  standard  point  kinetics  equations,  using  six  delayed  neutron 
groups  as  programmed  into  RELAP-3.  The  value  of  3 */i*,  the  effective 
delayed  neutron  fraction  over  pronyt  neutron  lifetime,  for  MH-1A  is  pro¬ 
vided.  Nonlinear  Doppler  feedback  in  the  form  of  reactivity  versus  fuel 
temperature  is  input  into  RE LAP- 3  along  with  constant  moderator  tenyerature 
and  void  reactivity  coefficients  corresponding  to  BOL  conditions.  The 
primary  shutdown  mechanism  during  the  LOCA  is  void  formation  in  the  core. 

For  conservatism  the  calculated  void  reactivity  coefficient  is  reduced  by 
25  percent  to  account  for  uncertainties .  The  spatial  distribution  of  the 
moderator  tenyerature  and  void  reactivity  feedback  is  assumed  to  be  pro¬ 
portional  to  the  product  of  the  square  of  the  axial  power  shape  and  the 
fraction  of  core  Vuluae  in  each  region. 

Heat  generated  during  blowdown  is  calculated  from  the  fission  and 
residual  heat- gene  rati  on  in  the  core  in  RELAP-3.  A  heat  balance  is  per¬ 
formed  in  each  ‘.xial  region.  Heat  transfer  rate  is  determined  based  on 
correlations  built  into  RELAP-3  and  calculated  local  conditions. 

6.  Safety  Injection.  A  provision  for  safety  injection  into  the 
reactor  \essel  downcomer  region  was  included  in  the  MH*1A  model.  A  constant 
flow  rate  of  100  gpui  was  assumed  consistent  with  the  design  value  for  the 
injection  system.  For  cold  leg  breaks,  it  was  assumed  that  ail  water 
injected  prior  to  the  end  of  blowdown  was  lost  out  of  the  break.  In  this 
context  the  end  of  blowdown  was  taken  as  the  time  when  zero  break  flow 

was  first  computed.  This  assunytion  has  negligible  effect  on  the  blowdown 
in  the  MH-1A  because  coolant  injection  on  the  order  of  100  gpm  is  insigni¬ 
ficant  compared  to  the  blowdown  flow  rates  for  large  breaks.  Details  of 
the  present  injection  system  design  are  presented  in  Appendix  A. 

C.  REFILL  CALCULATION 


1.  General.  The  refilling  phase  of  the  LOCA  is  begun  at  the  end  of 
blowdown  which  is  defined  as  the  time  when  zero  break  flow  is  first  computed. 
In  the  case  of  a  guillotine  break,  zero  break  flow  must  be  computed  from 
both  ends  of  the  break  before  the  end  of  blowdown  is  assumed.  After 
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receiving  an  actuation  signal,  refill  from  emergency  injection  tank  is 
initiated  when  the  reactor  Inlet  plenum  pressure  drops  below  240  psia. 

However,  in  accordance  with  the  AEC  interim  criteria,  all  coolant  injected 
prior  to  the  end  of  blowdown  is  assumed  lost  out  of  the  break.  After  the 
end  of  blowdown  the  coolant  is  assumed  tc  be  delivered  to  the  reactor 
vessel  inlet  plenum.  The  water  remaining  in  the  downcomer  region  below 
the  inlet  nozzle  and  in  the  lower  plenum  region  after  blowdown  is  con¬ 
servatively  neglected.  Refill  is  continued  until  the  coolant  reaches 
the  bottom  of  the  core  and  steam  is  generated  by  the  hot  fuel  rods  causing 
a  pressure  buildup  in  the  core  and  retarding  the  inlet  flooding  rate. 

Heat  transfer  is  evaluated  for  the  hot  fuel  rods  during  the  steam  cooling 
period  after  blowdown  when  the  core  is  uncovered.  After  core  reflooding 
is  initiated  heat  transfer  coefficients  from  the  FLSCRT  data  (Reference  12) 
can  be  used. 

2.  Vessel  Refill.  Driving  head  for  the  refill  injection  is  provided 
by  the  difference  between  the  nitrogen  pressure  in  the  tank  and  the  reactor 
vessel  back  pressure.  The  rate  of  coolant  injection  depends  on  the  magni¬ 
tude  of  the  driving  head  and  pressure  losses  in  the  piping  between  the 
injection  tank  and  reactor  vessel  downcomer.  The  gas  pressure  in  the 

tank  is  supplied  by  a  series  of  gas  bottles  and  is  regulated  such 
that  the  pressure  remains  constant  during  the  period  of  injection.  The 
nominal  design  injection  rate  Is  100  gpm  based  on  an  assumed  containment 
back  pressure  of  155  psia.  The  lower  containment  pressure  calculated  by 
RELAP-3  after  a  LOCA  should  result  in  a  higher  coolant  injection  rate. 

Figure  2  shows  the  refill  times  in  terms  of  the  emergency  injection  flow 
rate.  This  curve  is  based  on  the  reactor  vessel  volumes  and  a  constant 
coolant  injection  rate.  Boil  off  due  to  heat  transfer  from  the  reactor 
vessel  walls  and  hot  internals  is  not  included  in  the  calculations. 

An  estimate  of  the  maximum  temperature  rise  of  the  coolant  during 
refill  was  made  based  on  the  stored  energy  in  the  metal  below  the  core 
region  and  the  coolant  volume.  Assuming  complete  mixing  of  the  coolant, 
the  equilibrium  temperature  between  the  coolant  and  the  reactor  vessel  metal 
was  determined  to  be  185°F  with  an  initial  metal  temperature  of  500°F  and 
an  initial  water  temperature  of  85°F.  Therefore,  although  local  steam 
generation  would  occur  as  the  result  of  quenching  of  the  hot  metal  surfaces 
of  the  reactor  vessel,  the  bulk  coolant  temperature  would  remain  below  the 
saturation  temperature  (approximately  340° F)  while  the  coolant  level  is  below 
the  active  core  region.  For  rapid  refilling  rates,  a  relatively  small 
temperature  rise  in  the  bulk  coolant  would  be  expected.  This  is  desirable 
because  increased  subcooling  of  the  water  improves  the  flooding  heat  transfer. 

3.  Core  Reflooding.  Until  the  water  reaches  the  bottom  of  the  core, 
the  water  levels  in  the  downcomer  region  and  lower  plenum  beneath  the  core 
are  equal.  After  the  coolant  level  reaches  the  bottom  of  the  core  heat  is 
transferred  from  the  hot  fuel  rods  to  the  coolant.  Initially,  the  heat 
transferred  serves  to  remove  the  subcooling  of  the  water  entering  core 
and  the  core  flooding  rate  is  the  same  as  the  vessel  flooding  rate. 

However,  in  a  short  period  of  time  the  hot  fuel  rods  begin  to  generate 
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steam  In  the  core  region  which  must  be  vented  to  the  containment.  For 
the  hot  leg  breaks  this  steam  is  vented  directly  out  the  break  fr-.-m  the 
upper  plenum  region,  having  negligible  effect  on  the  core  flooding  rate. 
However,  for  the  cold  leg  break  the  steam  generated  in  the  core  region 
must  pass  through  the  entire  loop  before  being  vented  to  containment  out 
the  break.  The  flooding  rate  is  thus  limited  by  the  resistance  of  the 
loop  to  steam  flow. 

The  pressure  drop  along  the  flow  path  relieving  steam  after  the 
cold  leg  break  is  calculated  based  on  an  equivalent  loss  coefficient  for 
the  loop  corrected  for  existing  fluid  conditions.  The  reactor  coolant 
pump  rotors  are  assumed  to  be  locked  and  effects  of  superheating  in  the 
steam  generator  are  considered.  The  effect  of  steam  binding  tends  t< 
reduce  the  core  flooding  rate,  resulting  in  Increased  cladding  tempera¬ 
tures.  The  flooding  rate  during  the  steam  binding  period  is  a  function 
of  the  elevation  head  between  the  water  levels  in  the  downcomer  and  core 
regions,  the  loop  resistance,  and  the  fraction  of  the  inlet  flow  rate 
which  is  vaporized  and  must  be  relieved. 


For  the  purpose  of  calculating  the  relievable  steam  flow  rate, 
the  water  level  in  the  downcomer  is  assumed  to  have  reached  the  bottom 
of  the  inlet  nozzle  and  the  driving  force  for  the  steam  venting  is  assumed 
to  be  the  elevation  head  in  the  dowucomer  region  above  the  coolant  level 
in  the  core.  Table  4  shows  the  loss  coefficients  used  in  determining 
the  loop  equivalent  loss  coefficient.  The  relievable  steam  flow  rate 
defined  in  terms  of  the  equivalent  loss  coefficient  is  expressed  as  follows: 


Ws2  -  2r  Ps(P&  %) 
K/Az 


where :  W8 

Hd 

Pi 

Ps 

g 

A 

K 


relievable  steam  flow,  lb/sec 
downcomer  elevation  head,  ft 
downcomer  coolant  density,  lb/ft3 
saturated  steam  density,  lb/ft^, 
gravitational  constant,  ft/sec2 
primary  pipe  flow  area,  ft2 
equivalent  loss  coefficient 
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TABLE  4 


LOSS  COEFFICIENTS  FOR  MH-1A  COLD  LEG  STEAM  VENTING  CALCULATION 


Normal 

Operation 

Saturated 

Steam 

Superheated 

Steam 

Upper  Plenum  to  steam  generator  inlet 

1.62 

1.62 

1.62 

Steam  generator  inlet  to  break 

3.68 

4.67 

5.79 

Pumps  (locked) 

— 

6.10 

7.56 

Expansion  loss  to  containment 

— 

.98 

1.20 

TOTAL  equivalent  K 

5.30 

13.37 

16.17 

The  foilwing  assumptions  are  made  in  calculating  the  system 
loss  coefficient  and  relievable  steam  flow  rate: 

a.  The  reactor  coolant  pump  rotors  are  locked. 

b.  Steam  passing  through  the  steam  generator  is  superheated  to 

500°F. 

c.  The  coolant  leaving  the  core  is  saturated  steam. 

d.  The  containment  back  pressure  is  equal  to  the  initial  pre¬ 
break  pressure  plus  80  percent  of  the  increase  in  pressure  calculated  for 
the  accident . 

e.  The  emergency  injection  flow  does  not  affect  the  relievable 
steam  flow  or  pressure  drop  because  the  injection  does  not  occur  in  the 
steam  venting  flow  path  and  therefore,  no  flow  blockage  is  assumed. 

f.  The  additional  driving  pressure  for  core  reflooding  caused 
by  nitrogen  discharge  into  the  inlet  plenum  after  the  injection  tanks 
are  empty  is  conservatively  neglected. 

Figure  3  shows  the  relievable  steam  flow  rate  for  the  cold  leg 
breaks  at  the  pump  discharge.  The  actual  core  flooding  rate  would  depend 
upon  the  fraction  of  coolant  entering  the  core  which  is  vaporized  and  must 
be  relieved.  This  effect  should  be  calculated  based  on  the  nominal  flooding 
heat  transfer  data  from  Reference  12.  For  the  purpose  of  estimating  steam 
limited  core  flooding  rates  it  was  assumed  that  60  percent  of  the  coolant 
entering  the  core  was  vaporized  or  entrained.  Based  on  this  assumption,  the 
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steam  limited  flooding  rate  would  be  approximately  0.75  in/sec.  It  should 
be  pointed  out  that  a  significant  portion  oi  the  core  may  be  recovered 
before  steam  generation  causes  a  reduced  flooding  rate.  In  the  MH-1A 
the  hot  spot  at  BOL  occurs  at  approximately  8  inches  above  the  bottom  of 
the  core.  This  would  probably  be  effectively  cooled  before  steam  limiting 
could  occur.  Therefore,  the  most  important  factor  during  the  refill 
phase  appears  to  be  the  extent  of  core  cooling  during  the  period  of  time 
from  the  end  of  blowdown  to  the  time  the  injection  water  reaches  the  bottom 
of  the  core.  Heat  transfer  during  this  period  is  discussed  in  Section  VI-A 
with  respect  to  the  required  Emergency  Core  Cooling  System  performance. 

D.  FUEL  TEMPERATURE  TRANSIENT 

1.  General.  The  fuel  and  cladding  temperature  calculations  are  performed 
with  the  THETAl-3  computer  program  (Reference  8).  This  program  solves  the 
two  dimensional  transient  heat  conduction  equations  in  the  fuel  rod  and 

one  dimensional  conservation  of  energy  equations  for  the  associated  coolant 
channel.  The  surface  heat  transfer  model  in  THETAl-B  provides  the  boundary 
conditions  required  during  the  solution  of  the  fuel  rod  and  coolant  equations. 
The  required  input  Includes  coolant  channel  flow,  primary  system  pressure, 
core  power  generation  and  coolant  enthalpy  versus  time  as  obtained  from 
RELAP-3.  Following  blowdown,  the  heat  transfer  coefficient  versus  time 
during  reflooding  is  input.  Heat  generation  in  the  clad  due  to  metal  water 
reaction  is  also  included  in  the  model.  This  section  discusses  the  basic 
assumptions  and  methods  adopted  in  using  THETAl-B  for  the  analysis  of  the 
fuel  and  clad  temperature  transient. 

2.  Fuel  Red  Representation  The  fuel  rod  in  THETAl-B  is  represented 

by  six  equal  thickness  radial  regions  in  the  fuel,  two  regions  in  the  clad¬ 
ding  and  eleven  axial  regions  in  the  coolant.  Heat  transfer  in  the  axial 
direction  is  assumed  in  the  cladding,  but  is  conservatively  neglected  in 
the  fuel.  The  gas  gap  between  the  fuel  and  cladding  is  treated  as  a  heat 
conducting  path  with  negligible  heat  capacity.  The  gap  conductance  is 
calculated  based  upon  differential  expansion  of  the  fuel  and  cladding  and 
the  thermal  conductivity  of  the  gas.  The  initial  steady-state  fuel  tem¬ 
perature  distribution  in  THETAl-B  is  consistent  with  the  MH-1A  Core  4 
calculated  temperature  distribution  at  the  hot  spot  (Reference  13). 

3.  Heat  Transfer  Coefficients  During  Blowdown.  During  the  blowdown 
portion  of  a  LOCA  the  fuel  rod  surface  experiences  rapidly  changing  conditions 
of  flow,  pressure,  quality,  and  ^mperature.  The  heat  transfer  regime  and 
surface  heat  flux  are  determined  n  THETAl-B  on  the  basis  of  local  fluid 
conditions  and  cladding  surface  temperature.  The  heat  transfer  regimes 
which  could  be  encountered  during  a  LOCA,  and  the  corresponding  correla¬ 
tions  used  in  THETAl-B,  are  listed  in  Table  5.  Except  for  the  stable  film 
boiling  regime,  these  are  the  correlations  recommended  in  the  THETAl-B 
ceport.  For  film  boiling  heat  transfer  the  Groeneveld  correlation  (equa¬ 
tion  5.7  of  AECL-3281,  December  1969,  Reference  14)  was  used.  The  time  to 
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departure  from  nucleate  boiling  was  calculated  using  the  recommended 
critical  heat  flux  option  in  the  code  and  heat  transfer  after  departure 
from  nucleate  boiling  was  calculated  with  the  programmed  transition 
boiling  correlation.  Further  details  of  the  heat  transfer  model  in 
THETA1-B  are  discussed  in  Reference  8. 


TABLE  5 


HEAT  TRANSFER  CORRELATIONS 

DURING  BLOWDOWN  IN  THETA1-B 

Heat  Transfer 

Regime 

Correlation 

Reference 

1. 

Forced  convection  in 
sub cooled  liquid 

Dittus-Boelter 

15 

2. 

Nucleate  boiling 

Thom 

16 

3. 

Forced  convection 
vaporization 

S  chrock-Gros  sman 

17 

4. 

Transition  boiling 

McDonough,  Milich  and  King 

18 

5. 

Stable  film  boiling 

Groeneveld  (equation  5.7) 

14 

6. 

Pool  film  boiling 

Berenson 

29 

7. 

Forced  convection  in 
superheated  steam 

Dittus-Boelter 

15 

8. 

Critical  heat  flux 
correlations 

Combination  of  B&W-2 

Barnett  and  Modified  Barnett 

19,20,21 

4.  Hot  Channel  Input  Data.  The  boundary  conditions  used  in  solving 
the  fluid  energy  equations  in  THETA1-B  are  the  system  pressure,  flow  rate, 
and  enthalpy  of  the  coolant  entering  the  core  and  the  fuel  rod  surface  heat 
flux.  The  first  three  of  these  boundary  conditions  are  obtained  directly 
from  RELAP-3.  If  the  core  flow  at  the  upper  and  lower  core  junctions  in 
RELAP-3  is  in  the  same  direction  then  0.8  of  the  RELAP-3  smoothed  flow  at 
the  junction  entering  the  core  is  used.  If  the  flows  are  opposed,  a  stag¬ 
nation  point  in  the  core  is  assumed  and  zero  flow  is  used.  The  pressure 
and  enthalpy  used  in  THETA1-B  are  the  entering  plenum  conditions  obtained 
from  RELAP-3.  The  heat  flux  at  the  fuel  rod  surface  is  computed  by  the 
heat  transfer  routines. 
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5.  Reactor  Power.  The  reactor  power  in  the  fori?  of  internal  heal 
generation  as  a  function  of  time  is  input  to  THETAI-3.  The  normalized 
power  calculated  in  RELAP-3  from  the  kinetics  equations  considering  the 
various  reactivity  feedback  effects  is  used  until  shutdown  to  the  residual 
level  occurs.  Then  the  decay  heat  curve  described  in  the  proposed  ANS 
Standard  (Reference  22)  is  used.  This  curve  is  increased  by  20  percent 

to  account  for  uncertainties  and  includes  an  additional  allowance  for  heat 
generation  due  to  heavy  element  decay.  The  fraction  of  decay  heat  generated 
in  the  hot  fuel  rod  is  considered  to  be  100  percent  of  this  value. 

6.  Core  Power  Distribution.  The  MH-1A  Core  4  axial  power  distribution 
having  a  maximum  axial  peak  of  1.98  wa3  assumed  in  THETA1-B.  The  Core  4 
radial  peaking  factor  of  1.95  was  assumed  resulting  in  a  maximum  three 
dimensional  nuclear  peaking  factor  of  3.86.  An  additional  10  percent 
uncertainty  in  the  nuclear  power  distribution  was  assumed  and  a  3  percent 
uncertainty  due  to  fabrication  and  engineering  tolerances  was  added.  This 
results  in  a  peak  linear  heat  rate  of  19.1  kw/'ft  rated  power  and  21.1  kw/ft 
at  110  percent  power. 

7.  Fuel  Rod  Properties.  The  variation  of  conditions  in  the  core  during 
a  loss  coolant  accident  requires  a  knowledge  of  the  high  temperature  physical 
properties  of  the  core  materials.  Property  data  including  heat  capacity, 
thermal  conductivity,  density  and  mean  thermal  expansion  over  the  range  of 
temperatures  encountered  during  a  LOCA  are  incorporated  into  THETA1-B 
consistent  with  the  information  reported  in  Reference  23.  The  thermal 
conductivity  for  the  cladding  was  taken  as  th'-.  recommended  value  for 
stainless  steel  with  a  5  percent  uncertainty.  The  cladding  heat  capacity 
was  determined  from  the  weighted  contributions  of  the  principal  constituents 
of  the  stainless  steel  alloy.  The  property  data  for  the  UO2  fuel  used  in 
the  analysis  is  consistent  with  the  irradiated  UO2  data  reported  in  Reference 
23.  Further  details  of  the  materials  property  data  in  THETA1-B  are  presented 
in  Appendix  F. 

8.  Metal  Water  Reaction.  Transient  clad  surface  temperature  and  amount 
of  stainless  steel  reaction  is  determined  for  a  spectrum  of  peaking  factors 
from  the  THETA1-B  code.  The  stainless  steel-steam  reaction  rate  was  conserva¬ 
tively  taken  as  100  percent  of  the  value  determined  from  the  parabolic  rate 
law  (assuming  no  steam  limiting).  The  rate  equation  Incorporated  into 
THETA1-B  is  discussed  in  Appendix  G.  The  fraction  of  the  core  associated 
with  a  given  3  dimensional  peaking  factor  is  determined  by  combining  the 
output  of  the  two  dimensional  radial  peaking  factor  with  the  design  axial 
power  distribution  to  give  the  core  wide  power  peaking  distribution.  The 
fraction  of  the  core  associated  with  a  given  peaking  factor  is  then  combined 
with  the  fraction  of  metal  water  reaction  for  a  given  peaking  1  actor  to 
establish  the  percent  metal  water  reaction  in  the  core. 
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V.  RESULTS 


Analyses  were  performed  for  the  full  range  of  break  sizes  shown  in 
Table  1  using  the  conservative  assumptions  and  procedures  described  in 
Section  IV.  The  blowdown  portion  of  the  analysis  was  carried  out  cc-  che 
time  when  zero  break  flow  was  first  computed.  Information  regarding  the 
performance  of  the  ECCS  and  the  variation  of  significant  parameters  during 
blowdown  was  obtained.  Results  for  each  of  the  cases  are  sumarized  in 
the  form  of  curves  of  significant  parameters  as  a  function  of  time.  The 
fuel  temperature  transient  during  blowdown  and  during  the  initial  phase 
of  core  heat  up  was  calculated  for  the  hot  fuel  rod  using  THETA1-B.  Several 
representative  cases  are  discussed  in  detail  and  the  effect  of  continued 
core  heat  up  with  no  emergency  cooling  is  evaluated.  Results  of  this 
analysis  are  used  to  identify  the  most  severe  break  condition  and  determine 
the  minimum  time  for  core  cooling  to  become  effective. 

A.  BLOWDOWN 


The  sequence  of  events  occurring  during  the  subcooled  and  saturated 
phases  of  blowdown  of  the  primary  system  were  calculated  using  RELAP-3. 

The  results  of  these  analyses  for  the  spectrum  of  break  sizes  and  locations 
are  sumnarized  in  Figures  4  through  12.  These  curves  plot  the  variation 
of  inportant  parameters  during  the  depressurization  including  the  break 
flow  rate,  the  system  pressure,  reactor  power,  core  flow  rate,  and  inlet 
enthalpy  as  a  function  of  l  me.  For  the  cold  leg  breaks  in  which  the  flow 
rate  through  the  core  reverses,  the  upper  plenum  conditions  are  presented. 

As  shown  in  the  figures,  the  system  is  rapidly  depressurized  immediatelv 
after  the  break  until  the  vapor  pressure  of  the  coolant  is  reached.  During 
this  subcooled  phase  pressure  waves  are  propagated  through  the  system  at 
the  speed  of  sound  resulting  in  significant  pressure  differentials  across 
system  coaponents.  Figure  13  and  14  show  the  calculated  core  pressure  drop 
as  a  function  of  time  for  the  breaks  at  the  inlet  and  outlet  of  the  reactor 
vessel  indicating  that  pressure  differentials  between  80  and  90  psi  can 
occur  momentarily  across  the  core.  However,  the  pressure  oscillatioi  are 
quickly  daoped  due  to  the  formation  of  voids  in  the  system.  Within 
approximately  0.02  seconds  the  initial  subcooled  depressurization  is  com¬ 
plete  and  saturated  blowdown  conmences  with  choked  flow  occurring  at  the 
break  location.  The  core  power  level  is  rapidly  reduced  due  to  formation 
of  voids  in  the  core  region  and  the  flow  rate  stabilizes  and  decreases. 
Reactor  scram  due  to  low  pressure  occurs  after  approximately  two  seconds, 
being  delayed  by  the  blowdown  of  the  pressurizer  liquid  volume  through  the 
surge  line.  Heat  transfer  in  the  core  during  blowdown  continues  at  a  high 
rate  until  the  critical  heat  flux  is  reached  at  which  time  a  significant 
reduction  in  heat  transfer  at  the  clad  surface  occurs.  For  large  breaks  in 
the  HH-1A  the  period  of  effective  core  cooling  is  from  three  to  five  seconds 
and  the  blowdown  is  essentially  complete  alter  approximately  si\  set  ads. 
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Table  6  shows  the  tides  at  which  significant  events  occur  during  the 
blowdown  as  calculated  by  RE LAP-3.  The  table  includes  the  tiae  at  which 
zero  break  flow  is  first  co^mted,  the  tines  at  which  low  pressurizer  pres¬ 
sure  and  high  containment  pressure  actuation  signals  are  initiated,  and 
the  time  at  which  the  pressure  in  the  inlet  plenum  drops  below  the  pressure 
in  the  safety  injection  tank  permitting  the  flow  of  emergency  coolant  into 
the  vessel.  As  shown  in  Table  6,  the  delay  in  the  receipt  of  the  safety 
injection  actuation  signal  from  the  pressurizer  is  the  controlling  signal 
for  safety  injection  for  large  breaks.  The  high  containment  pressure  set- 
point  is  reached  shortly  after  the  breaks  occurs.  Actuation  of  the  valves 
in  the  nitrogen  pressurizing  and  safety  injection  lines  and  flow  of  coolant 
in  the  piping  from  the  check  valve  to  the  injection  point  would  be  a  source 
of  additional  delay. 

The  values  of  various  calculated  parameters  at  the  time  when  zero  break 
flow  first  occurs  are  shown  in  Table  7.  The  results  of  the  calculation 
with  RELAP-3  as  shown  indicate  that  from  96  to  98  percent  of  the  initial 
primary  system  coolant  is  lost  out  of  the  break  by  the  end  of  blowdown  and 
the  liquid  mass  remaining  in  the  reactor  inlet  plenum  below  recirculation 
nozzles  is  a  small  percentage  of  the  initial  inlet  plenum  volume.  These 
results  ate  consistent  with  experimental  observations  of  coolant  remaining 
in  the  system  after  rupture  and  decompression  of  the  semiscale  test  loop 
(Reference  30).  The  calculated  containment  pressure  after  blowdown  varies 
from  approximately  118  to  128  psia  considerably  below  the  design  pressure 
of  155  psia.  Smaller  break  sizes  resulting  in  a  larger  amount  of  energy 
transferred  from  the  core  to  the  containment  could  result  in  higher  con¬ 
tainment  back  pressure  after  blowdown.  However,  for  the  break  sizes 
analyzed  the  lower  containment  back  pressure  would  indicate  an  improvement 
in  safety  Injection  flow  rate  of  approximately  10  percent  over  the  nominal 
design  value.  Also  shown  is  the  core  average  fuel  tenperature  in  the  high 
heat  flux  axial  region  at  the  end  of  blowdown.  These  temperatures  indicate 
the  relative  heat  generation  and  removal  rates  for  the  average  core  during 
blowdown  and  depend  primarily  upon  the  length  of  time  of  blowdown  and  the 
effective  heat  transfer  coefficient  calculated  in  RELAP-3.  These  results 
appear  to  confirm  that  the  double-ended  guillotine  break  of  the  primary  pipe 
is  the  most  severe  break  with  respect  to  the  clad  temperature  transient 
because  of  lower  amount  of  energy  removed  during  blowdown. 

The  effect  of  the  reactor  coolant  punp  characteristics  during  blowdown 
was  evaluated  for  the  0.6  DE  slot  break  in  the  cold  leg.  Two  rases  were 
examined;  the  first  considered  a  break  at  the  pump  discharge  with  assumptions 
resulting  in  minimum  pump  resistance  and  the  second  case  considered  a  break 
at  the  pump  inlet.  For  the  break  at  the  pump  discharge  it  was  assumed 
that  the  pumps  remained  running  during  the  blowdown,  following  the  pump  head 
curve  until  cavitation  occurred.  For  the  break  at  the  pump  inlet  it  was 
assumed  that  the  pumps  tripped  at  the  time  the  break  occurred  due  to  a  loss 
of  power  and  that  a  resistance  due  to  a  free  spinning  rotor  existed  after 
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0.5  ft2  Break  at  Pump  Discharge  (S)  8.70  0.41  5.18  6.35 

.18  ft2  Surge  Line  Break  19.18  0.52  3.13  16.13 


CALCULATED  PARAMETERS  AT  THE  END  OF  BLOWDOWN  (BREAK  FLOW 
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the  pumps  shutdown.  Table  8  summarizes  the  assumptions  used  in  the  two 
cases.  Figures  15  and  16  show  the  effects  of  the  assumptions  on  the  re¬ 
sistance  across  the  punp  and  the  flow  rate  through  the  break.  The  more 
rapid  blowdown  for  the  case  with  the  minimum  pump  resistance  (case  1) 
resulted  in  reduced  heat  removal  from  the  core  and  consequently  the  more 
severe  temperature  transient.  These  assumptions  were  used  for  the  re¬ 
mainder  of  the  cold  leg  breaks . 


TABLE  8 

REACTOR  COOLANT  PUMP  ASSUMPTIONS 


Case 

Break  Location 

Pump  Trip 

Pump  Head 
after 

Cavitation 

Shutdown 

Pump 

Resistance 

1 

Pump  Discharge 

No 

None 

None 

2 

Pump  Inlet 

Yes 

None 

Free  Rotor 

B.  CORE  HEAT  UP  WITH  PRESENT  INJECTION  SYSTEM 


The  thermal  response  of  the  fuel  and  the  cladding  at  the  hot  spot  is 
calculated  during  blowdown  and  Initial  core  heat  up  using  THETA1-B. 

Initially  after  the  break  occurs,  high  convective  and  nucleate  boiling 
heat  transfer  rates  continue  to  prevail  in  the  core.  However,  as  the 
blowdown  proceeds  conditions  resulting  in  departure  from  nucleate  boiling 
occur  and  the  surface  heat  flux  is  drastically  reduced.  After  the  blow¬ 
down  is  complete,  the  majority  of  the  coolant  inventory  is  lost  and  primary 
mechanism  for  removing  heat  from  the  core  is  through  radiation  and  convec¬ 
tion  to  steam  generated  from  coolant  in  contact  with  vessel  walls  beneath 
the  core. 

The  present  MH-1A  ECCS  is  designed  to  deliver  coolant  at  the  rate  of 
100  gpm,  resulting  in  a  time  of  approximately  6  minutes  to  refill  the  reactor 
vessel  to  the  bottom  of  the  core.  Thus,  with  the  present  system  there 
exists  an  extended  period  of  time  in  which  steam  cooling  is  relied  upon. 

For  the  purpose  of  estimating  the  minimum  time  required  before  emergency 
cooling  must  become  effective  and  of  assessing  the  relative  severity  of  the 
various  cases,  the  assumption  of  no  heat  transfer  at  the  cladding  surface 
was  applied  at  the  end  of  blowdown.  It  is  clear  that  this  assumption  is 
conservative  and  results  in  continued  core  heat  up  at  a  rate  depending 
primarily  on  the  stored  energy  in  the  fuel  at  the  end  of  blowdown  and 
the  residual  heat  generation  rate.  At  temperatures  above  approximately 
2100°F,  additional  heat  is  added  due  to  the  reaction  of  cladding  with  steam. 
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Figures  17  through  19  show  the  results  of  the  hot  spot  temperature 
transient  during  blowdown  and  the  initial  core  heat  up  for  the  hot  and 
cold  leg  double  ended  guillotine  breaks  and  the  pressurizer  surge  line 
break.  The  boundary  conditions  from  RELAP-3  were  used  up  until  the  end 
of  blowdown  at  which  time  zero  heat  transfer  was  assumed  at  the  clad 
surface.  Fuel  and  clad  temperatures,  heat  transfer  coefficient,  fluid 
temperature,  and  local  quality  and  void  fraction  at  the  location  of  the 
hot  spot  are  plotted  for  each  case.  The  plots  of  fuel  and  clad  tenperature 
variation  with  time  illustrate  the  redistribution  of  sensible  heat  in  the 
fuel  and  cladding  immediately  following  the  rupture.  This  effect  results 
in  the  major  source  of  clad  heating  up  to  about  20  seconds  after  the  break. 
Following  this  time  continued  core  heat  up  due  to  decay  heat  generation 
and  metal  water  reaction  occurs. 

The  effect  of  continued  cladding  heat  up  following  the  double-ended 
hot  leg  break  for  specified  volume  percentages  of  the  core  is  shown  in 
Figure  20.  The  times  f~r  the  cladding  to  reach  2300°F  and  the  melting 
temperature  of  stainlet  steel  (2550°F)  are  shewn  in  Table  9.  The  calcu¬ 
lation  performed  for  th.s  curve  assumed  that  there  was  no  emergency  cooling 
after  blowdown  but  that  sufficient  steam  was  available  for  metal  water 
reaction.  Redistribution  of  the  initially  peaked  temperature  profile  occurs 
in  approximately  20  seconds  after  which  continued  core  heat  up  occurs  at  a 
rate  depending  on  the  decay  heat  level  and  local  power  peaking.  Figure  21 
plots  the  percent  of  metal  water  reaction  as  a  function  of  time  for  the 
hot  spot  and  overall  core.  From  this  figure  it  can  be  seen  that  the  metal 
water  reaction  becomes  significant  at  temperatures  above  approximately 
2100°F  and  rapidly  increases  as  the  temperature  rises. 

Table  10  presents  the  results  of  the  clad  temperature  calculation  for 
each  of  the  break  conditions.  The  values  shown  are  for  30  seconds  after 
rupture  at  which  time  redistribution  of  the  clad  and  fuel  temperature  is 
essentially  complete.  The  calculations  conservatively  assumed  no  emergency 
cooling  after  blowdown.  From  these  results  the  cold  leg  break  appears  to 
be  in  general  more  severe  than  the  hot  leg  break  due  to  the  lower  heat 
transfer  during  blowdown.  The  amount  of  heat  removed  from  the  hot  spot 
during  blowdown  largely  determines  the  period  of  time  remaining  before  the 
emergency  cooling  must  become  effective.  On  this  basis,  the  minimum  delay 
for  the  ECCS  to  become  effective  is  approximately  30  seconds  after  the 
break.  This  delay  results  from  the  conservative  assumption  that  unlimited 
steam  would  be  available  for  the  metal  water  reaction  but  not  for  cooling 
the  core.  In  reality  the  actual  delay  time  would  be  longer  depending  upon 
the  amount  of  heat  transfer  occuring  due  to  radiation  and  convection  after 
blowdown.  The  effect  of  heat  transfer  at  the  hot  spot  during  this  period 
is  discussed  in  Section  VI-A. 
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TABLE  9 


TIME  REQUIRED  FOR  CLADDING  TO  REACH  2300°F  AND  MELTING  TEMPERATURE 
OF  STAINLESS  STEEL  FOR  DE  HOT  LEG  BREAK  WITH  NO  EMERGENCY  COOLING 


Percent  of  Time  to  Reach  Time  to  Reach  Melting 

Clad  Volume  2300°F,  sec  Temperature,  sec 


0  (hot  spot) 

44 

65 

2 

78 

101 

10 

129 

155 

25 

250 

300 

TABLE 

10 

RESULTS  OF  CLAD  TEMPERATURE 

BREAK 

TRANSIENT  AT  HOT  SPOT 
NO  EMERGENCY  COOLING 

30  SECONDS  AFTER 

Break  Type 

Hot  Spot 

Clad  Temp 
°F 

%  Metal 

Water  Reaction 
(Local) 

DE  Hot  Leg  Break  (G) 

2091 

.56 

DE  Hot  Leg  Break  (S) 

2075 

.49 

SE  Hot  Leg  Break  (S) 

1796 

.02 

DE  Break  at  Pump  Discharge  (G) 

2290 

2.27 

DE  Break  at  Pump  Discharge  (S) 

2163 

1.00 

.6  DE  Break  at  Pump  Discharge  (S) 

2275 

2.08 

.6  DE  Break  at  Pump  Inlet  (S) 

2167 

.98 

0.5  ft2  Break  at  Punp  Discharge 

(S) 

2284 

2.17 

0.18  ft2  Surge  Line  Break 

969 

0 
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VI.  REQUIRED  EMERGENCY  CORE  COOLING  SYSTEM  PERFORMANCE 


This  section  discusses  the  design  and  performance  of  the  MH-1A  ECCS 
on  the  basis  of  the  calculated  results  of  the  LOCA  and  AEC  Interim  Acceptance 
Criteria,  and  estimates  the  required  cooling  capability  of  the  system. 
Consideration  of  the  post  blowdown  environment  in  the  vessel  and  heat  sources 
in  the  core  pertaining  to  the  effectiveness  of  the  ECCS  are  outlined. 

The  primary  requirement  for  the  MH-1A  ECCS  is  to  provide  reliable  and 
effective  core  cooling  within  approximately  30  seconds  after  a  large  rupture 
in  the  primary  piping  by  rapidly  adding  coolant  into  the  vessel  and  core 
region.  For  the  MH-1A  the  loss  of  coolant  inventory  after  a  large  break 
occurs  within  approximately  four  to  six  seconds.  This  rapid  blowdown  is 
due  primarily  to  the  potentially  large  break  size  in  a  single  loop  system 
relative  to  the  vessel  size  and  coolant  inventory  in  the  system.  However, 
in  a  single  loop  system  a  large  fraction  of  the  coolant  is  forced  to  pass 
through  the  core  region  during  blowdown  and  relatively  large  heat  transfer 
rates  exist  in  the  core.  After  blowdown  the  ECCS  must  be  capable  of  con¬ 
trolling  the  decay  heat  generation  until  the  cladding  is  quenched  thereby 
reducing  its  temperature  to  that  of  the  coolant.  The  semi-open  lattice  design 
of  the  MH-1A  should  promote  good  distribution  of  the  coolant  through  the 
core  region  before  and  during  core  flooding.  After  the  core  is  quenched 
long  term  core  cooling  would  be  provided  by  the  decay  heat  removal  system. 

The  primary  requirement  for  this  system  is  reliability  due  to  the  extended 
period  of  time  it  would  be  expected  to  operate  following  a  LOCA. 

A.  HEAT  REMOVAL  CAPABILITY 


Following  a  rupture  in  the  primary  system  the  ECCS  must  function  such 
that  the  fuel  rod  heat  energy  is  removed,  and  uncontrolled  cladding  heat  up 
is  prevented.  The  major  sources  of  energy  during  the  LOCA  include  the 
initial  stored  energy  in  the  fuel,  the  fission  energy  generated  prior  to 
complete  shutdown  from  voids  or  control  rod  motion,  the  residual  heat 
generation  due  to  fission  product  decay,  and  possible  metal  water  reaction 
energy. 

The  heat  transfer  during  blowdcwn  removes  a  portion  of  the  stored  and 
fission  energy  in  the  fuel,  thereby  lowering  the  average  fuel  temperature 
at  the  hot  spot.  Following  a  lupture  in  the  primary  system,  redistribution 
of  sensible  heat  occurs  such  that  the  cladding  and  fuel  temperatures  approach 
a  common  temperature  approximately  equal  to  the  fuel  average  temperature 
after  blowdown.  For  the  MH-1A  the  blowdown  removes  approximately  600° F 
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average  tenperature  from  the  hot  spot.  Following  blowdown  ti.s.  major  source 
of  heat  in  the  fuel  is  due  to  fission  product  decay.  The  heat  flux  which 
must  be  removed  from  the  hot  spot  and  average  core  as  a  function  of  time 
after  shutdown  is  shown  in  Figure  22.  This  curve  is  based  upon  the  ANS 
Standard  (Reference  22)  for  decay  heat  and  the  MH-1A  power  peaking  factors. 
Above  temperatures  of  approximately  2100°F  the  metal  water  reaction  begins 
to  become  significant  i.  Iding  additional  heat  which  must  be  removed.  On 
this  basis  the  heat  transfer  coefficient  necessary  to  prevent  continued 
heat  up  of  the  fuel  and  cladding  above  2100°F,  assuming  a  sink  temperature 
corresponding  to  saturated  steam  (340°F) ,  would  be  approximately  15  Btu/hr 
ft2°F. 

Figures  23  and  24  show  the  hot  spot  cladding  temperature  transients 
for  the  double  ended  hot  and  cold  leg  guillotine  breaks  assuming  variable 
heat  transfer  coefficient  at  the  end  of  blowdown  and  a  sink  temperature 
of  340°F.  The  actual  heat  transfer  at  the  end  of  blowdownwill  depend  on 
the  calculated  local  conditions  in  the  core  region  and  can  be  determined 
using  correlations  for  forced  convection  to  steam  (References  26  and  27) 
and  radiation  to  steam  and  adjacent  fuel  rods  (Reference  28).  Other 
considerations  would  be  the  amount  of  steam  available  in  the  core  region, 
cladding  surface  temperatures,  and  type  of  fluid  flow  existing  in  the 
coolant  channels.  Conservative  estimates  of  the  radiation  to  satu¬ 
rated  steam  result  in  heat  transfer  coefficients  of  approximately  3  to  5 
Btu/hr  ft2°F  and  the  minimum  convective  heat  transfer  coefficient  due  to 
laminar  flow  is  approximately  3.5  Btu/hr  ft2oF.  Thus,  if  sufficient  steam 
were  provided  in  the  core  region  following  blowdown,  heat  transfer  coefficient 
of  6  to  9  or  larger  could  be  expected.  This  would  delay  the  time  before 
effective  cooling  of  the  core  is  required  by  an  additional  20  to  40  seconds 
or  longer. 

After  core  flooding  is  initiated  heat  transfer  from  the  FLECHT  data 
can  be  applied  (Reference  12).  The  major  differences  between  the  FLECHT 
tests  and  the  MH-1A  is  the  heated  length  of  the  fuel  rods  (12  feet  as  opposed 
to  3  feet  in  the  MH-1A)  and  the  elevation  of  the  elevation  of  the  peak  power 
density  (6  feet  in  the  FLECHT  tests  as  opposed  to  less  than  one  foot  in 
the  MH-1A) .  As  a  result  of  these  differences,  effective  cooling  of  the 
hot  spot  in  the  MH-1A  could  be  expected  shortly  after  the  coolant  reached 
the  bottom  of  the  core  due  to  frothing  and  splashing  of  the  coolant  caused 
by  steam  generation  in  the  flooded  region.  Detailed  calculations  of  the 
refill  and  reflooding  transient  considering  the  improved  system  design  will 
be  required  to  determine  the  actual  temperature  time  history  of  the  cladding 
and  the  performance  of  the  system  relative  to  the  AEC  Interim  Acceptance 
Criteria. 
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B.  DESIGN  CONSIDERATIONS 


In  general  the  design  Improvements  for  the  MH-1A  ECCS  should  be  con¬ 
cerned  with  improving  the  response  time  of  the  injection  system  and 
providing  sufficient  redundancy  and  reliability  in  the  system  to  assure 
effectiveness.  Specifically,  the  evaluation  of  the  MH-1A  ECCS  should 
consider  the  means  of  cooling  the  core  over  the  spectrum  of  breaks  from 
a  small  leak  to  doiile  ended  rupture  of  the  primary  pipe.  The  emergency 
coolant  for  large  breaks  should  be  added  in  such  a  way  that  effective 
cooling  of  the  core  region  could  be  accomplished  rapidly  after  the  break 
occurs.  Bottom  flooding,  as  exists  with  the  present  system,  and  flooding 
frop  above  the  core  could  be  considered.  Redundant  systems  would  be 
desirable  from  the  standpoint  of  the  single  failure  criterion.  Testing  of 
the  Injection  system  to  demonstrate  the  actual  flow  rates  achievable  would 
aid  in  the  evaluation  of  the  Improved  system.  The  effect  of  loss  of  power 
and  startup  delays  of  the  emergency  power  equipment  should  be  evaluated  and 
a  means  of  periodically  testing  the  operation  of  majo/  conponents  to 
assure  reliability  and  effectiveness  of  the  system  should  be  included. 
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FIGURE  1  .  MH-1A.  MODEL  -  SINGLE  ENDED  HOT  LEG  SLOT  BREAK 
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FIGURE  2  REFILL  TIMES  VERSUS  EMERGENCY  INJECTION  FLOW  RATE 


31 


9  November  1972 


FIGURE  4A  DOUBLE  ENDED  HOT  LEG  BREAK  (GUILLOTINE)  BREAK  FLOW  VS.  TIME 


FIGURE  4B  DOUBLE  ENDED  HOT  LEG  BREAK  (GUILLOTINE}  7RESSURE  VS.  TIME 


FIGURE  4C  DOUBLE  ENDED  HOT  LEG  BREAK  (GUILLOTINE)  POWER  LEVEL  VS.  TIME 
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FIGURE  5A.  DOUBLE  ENDED  HOT  LEG  BREAK  (SLOT)  BREAK  FLOW  RATE  VS.  TIME 


FIGURE  5B.  DOUBLE  ENDED  HOT  LEG  BREAK  (SLOT)  PRESSURE  VS 


'IGURE  5D .  DOUBLE  ENDED  HOT  LEG  BREAK  (SLOT)  CORE  ELOW  VS.  TIME 
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FIGURE  6B.  SINGLE  ENDED  HOI  LEG  BREAK  (SLOT)  PRESSURE  VS.  TIME 


TIME. SECONDS 

FIGURE  6C .  SINGLE  ENDED  HOT  LEG  BREAK  (SLOT)  POWER  LEVEL  VS.  TIME 
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FIGURE  6D.  SINGLE  ENDED  HOT  LEG  BREAK  (SLOT)  CORE  FLOW  VS.  TIME 
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FIGURE  7A.  DOUBLE  ENDED  BREAK  AT  pirMT>  nTRCHARCR  (GUILLOTINE)  BREAK 
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1.00  2.00  3.00  4.00  5.00  6-00  7.00 

TIME  .SECONDS 

FIGURE  7B .  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (GUILLOTINE) 

PRESSURE  VS.  TIME 


TGURE  7C.  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (GUILLOTINE)  POWER  LEVEL  VS.  TIME 


FIGURE  7D.  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (GUILLOTINE)  CORE  FLOW  VS.  TIME 
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FIGURE  «A.  DOUBLE  ENDED  COLD  LEG  AT  TUMP  DISCHARGE  (SLOT)  BREAK 
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FIGURE  8B.  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  PRESSURE  VS.  TIME 


FIGURE  8C.  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (ST GIT)  POWER  LEVEL  VS.  TIME 
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FIGURE  8D.  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  CORE  FLOW  VS.  TIME 


DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  INLET  ENTHALPY  VS.  TIME 
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DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  BREAK  FLOW  RATE  VS.  TIME 


FIGURE  9B.  DOUBLE  EN’DED  COLD  L'JG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  PRESSURE  VS.  TIME 
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T I  ME .SECONOS 

FIGURE  9D.  0.6  DOUBLE  ENDED  COLD  LEG  BREAK  AT  I’UMT  DISCHARG 


62  9  November  1972 


I 


TIME  .SECONDS 

FIGURE  1QA.  0.6  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  INLET  (SLOT)  BREAK  FLOW  RATE  VS.  TIME 


FIGURE  IOC.  0.6  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  INLET  (SLOT)  POWER  LEVEL  VS.  TIKE 
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TIME  . SECONDS 

FIGURE  lOF..  0.6  DOUBLE  ENDED  COLD  LEG  BREAK  AT  PUMP  INLET  (SLOT)  INLET  ENTHALPY  VS.  TIME 
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FIGURE  11A.  0.5  ft2  BREAK  AT  PUMP  DISCHARGE  (SLOT)  BREAK  FLOW  VS.  TIME 


rIGURE  11B .  0.5  ft2  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  PRESSURE 
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-EG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  POWER 
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FIGURE  HE.  0.5  ft2  COLD  LEG  BREAK  AT  PUMP  DISCHARGE  (SLOT)  INLET  ENTHALPY  VS.  TIME 
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FIGURE.  12  A.  0.18  ft2  PRESSURIZER  SURGE  LINE  BREAK  BREAK  FLOW  RATE  VS. 
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FIGURE  12B.  0.18  ft2  PRESSURIZER  SURGE  LINE  BREAK  PRESSURE  VS.  TIME 


FIGURE  12C.  0.18  ft2  PRESSURIZER  SURGE  LINE  BREAK  POWER  LEVEL  VS. 


PRESSURIZER  SURGE  LINE  BREAK  INLET  ENTHALPY  VS.  TIME 


FIGURE  14.  DOUBLE  ENDED  COLD  LEG  BREAK  (SLOT)  CORE  PRESSURE  DROP  VS.  TIME 
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FICURE  16.  0.6  DOUBLE  ENDED  COLD  LEG  BREAK  (SLOT)  BREAK  FLOW  KATE  VS.  TIME 
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FIGURE  17A  DOUBLE  ENDED  HOT  LEG  BREAK  (GUILLOTINE) 
CLAD  AND  FUEL  TEMPERATURE  AT  HOT  SPOT  V 
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FIGURE  17P  DOUBLE  ENDED  HOT  LEG  RRE 
HEAT  TRANSFER  COEF  FICIE.N" 


FIGURE  17D  DOUBLE  ENDED  HOT  LEG  BREAK  (GUILLOTINE) 

QUALITY  AND  VOID  FRACTION  AT  HOT  SPOT  VS.  TIME 
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FIGURE  18B  DOUBLE  ENDED  COLD  LEU  BREAK  ‘(GUILLOTINE) 

HEAT  TRANSFER  COEFFICIENT  AT  HOT  SPOT  VS.  TIME 
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TIME. SECONDS 

FIGURE  19B  0.18  ft2  SURGE  LINE  BREAK 

HEAT  TRANSFER  COEFFICIENT 
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IE  .SECONDS 

.8  ft2  SURGE  LINE  BREAK 

UlITY  AN'-)  VOIP  FRACTION  AT  HOT  SPOT 
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FIGURE  20.  MH-1A  PPl'PT.F  FVPFP  POT  LEG  BREAK  GUILLOTINE 
CLAD  TEMPERATURE  VS.  TIME 
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FIGURE  21.  METAL  WATER  REACTION  VS.  TIME 
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APPENDIX  A 

EMERGENCY  CORE  COOLING  SYSTEM 

The  Bsergency  Core  Cooling  System  presently  installed  in  the  MH-1A 
consists  of  emergency  injection,  coolant  charging  and  decay  heat  removal 
functions  performed  by  the  Coolant  Charging  and  Chemical  Addition  System. 
This  system  is  activated  in  the  event  of  a  loss  of  coolant  accident  so  as 
to: 


a.  Inject  borated  water  into  the  reactor  vessel  in  sufficient 
qiantlty  to  recover  the  core  and  to  render  the  reactor  subcritical  with 
all  rods  withdrawn. 


b.  Provide  emergency  decay  heat  removal  by  injecting  water  into  the 
reactor  vessel  immediately  following  loss  of  coolant. 

c.  Provide  makeup  with  coolant  charging  puaps  and  long  term  core 
cooling. 

A  brief  description  of  the  major  components  and  operation  of  the 
system  in  the  event  of  a  LOCA  is  presented  below: 

Emergency  Coolant  Injection.  The  emergency  injection  system  is  illustrated 
in  Figure  A-l.  This  system  consists  of  an  injection  tank,  twelve  bottles 
of  compressed  gas,  four  solenoid  valves,  two  check  valves  and  associated 
piping.  The  emergency  injection  tank  is  a  350  cubic  foot  tank  located 
on  the  platform  deck  in  the  reactor  access  compartment.  The  design  pres¬ 
sure  of  the  tank  is  235  psig  and  design  temperature  is  125°F.  The  injection 
tank  contains  a  minimum  of  1000  gallons  of  borated  water  which  is  sufficient 
coolant  volume  to  recover  the  core.  The  compressed  gas  is  contained  in 
twelve  bottles  located  on  the  upper  deck  in  the  gas  bottle  storage  room. 

The  bottles  are  arranged  in  two  banks  of  six,  each  bank  connected  to  the 
injection  tank  by  a  line  containing  a  solenoid  valve.  Upon  coincident 
signals  of  low  pressurizer  pressure  and  high  containment  pressure,  both 
solenoid  valves  open  allowing  the  gas  to  pass  through  the  pressure  regu¬ 
lating  valves  and  into  the  injection  tank,  thus  providing  the  necessary 
pressure  to  force  the  borated  water  from  the  tank  into  the  reactor  vessel. 
Redundant  solenoid  valves  assure  that  emergency  coolant  is  provided  if  one 
of  the  valves  in  the  pressurizing  or  coolant  injection  lines  should  fail. 

Coolant  Charging  and  Decay  Heat  Removal.  The  emergency  coolant  charging 
system  consists  of  a  10,000  gallon  emergency  water  storage  tank  located 
external  to  the  containment  vessel,  piping  and  valves  which  connect  the 
tank  to  the  reactor  vessel  fill  nozzle  line,  and  two  positive  displacement 
charging  pumps.  Upon  a  LOCA  actuation  signal  one  charging  pump  is  started 
automatically  from  the  450  volt  emergency  bus  providing  8.5  gpm  continuous 
makeup  flow  to  the  reactor  vessel  in  addition  to  emergency  coolant  flow  from 
the  injection  tank.  There  is  sufficient  water  in  the  emergency  water 
storage  tank  to  remove  the  equivalent  of  the  first  five  days  of  decay  heat 
generation.  Since  the  tank  water  can  be  replenished,  the  core  can  be  kept 
submerged  indefinitely.  Electrical  power  for  charging  pumps  is  backed  up 
by  three  independent  diesel  generators  and  the  emergency  diesel  as  well  as 
onshore  power  when  MH-1A  is  connected  to  a  dependable  power  grid. 
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PIGURE  A-l 

EMERGENCY  INJECTION  COOUNG  WATER  SYSTEM 
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APPENDIX  B 


PRIMARY  SYSTEM  GEOMETRY  AND  INITIALIZATION 


RELAP-3  enploys  a  one  dimensional  representation  of  the  primary'  system 
In  terms  of  fluid  volume  elements  and  junctions  describing  the  plenums 
piping,  reactor  core,  steam  generator,  pressurizer,  and  containment. 
Junctions  can  model  pumps  or  check  valves  and  can  simulate  leak  or  fill 
vater  sources.  Input  data  describes  the  interconnection  of  system  volumes, 
physical  dimensions  of  the  elements,  and  initial  steady  state  conditions. 
Table  B~1  briefly  summarizes  data  for  the  MH-1A  model  representing  a  single 
ended  hot  leg  riot  break  shown  in  Figure  1  of  the  main  text.  Other  models 
are  similar  except  for  break  type  and  location. 

Appropriate  data  describing  the  Initial  steady  state  conditions  in  the 
system  is  supplied  to  the  code.  The  initial  pressure  and  tenperatures 
distributions,  mass  flow  rate,  and  power  level  are  specified  assuming  that 
the  power  produced  in  the  core  is  being  removed  by  steam  generator. 

Friction  coefficient  data  required  for  the  solution  of  the  momentum  equation 
is  calculated  internally  in  RELAP-3  based  on  the  specified  input  pressure 
distribution  and  flow.  If  the  initial  flow  through  a  junction  is  zero,  then 
the  friction  coefficient  is  supplied.  Similarly  the  resistance  of  the  pumps 
after  shutdown  is  provided  as  discussed  in  Appendix  C. 

The  time  dependent  solution  is  determine.!  in  RELAP-3  by  integrating  the 
equations  describing  mass,  energy,  and  momentum  balances  for  each  of  the 
elements  using  the  thermodynamic  equations  of  state  to  obtain  the  fluid 
properties.  The  contents  of  volume  elements  during  blowdown  are  assumed  to 
be  homogeneous.  Fluid  properties  are  considered  uniform  with  the  pressure 
occurring  at  the  center  of  gravity.  Various  effects  such  as  pump  charac¬ 
teristics,  reactivity  feedback,  choked  flow,  core  and  steam  generator  heat 
transfer,  and  safety  injection  are  taken  into  account  during  the  transient 
solution . 
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TABLE  B-l 


SUMMARY  OF  RELAP-3  INPUT  DATA  FOR  MH-1A  MODEL 
SINGLE  ENDED  HOT  LEG  BREAK 


Volume  Data 


Volume  Volume  Elevation*  Initial 

No.  Descriptio.i  (ft3)  (ft)  Pressure 


(psia) 

1 

Reactor  Inlet  Plenum 

64 

0.0 

1508 

2 

Core  Volume  1 

7,3 

4.41 

1507 

3 

Core. Volume  2 

7.3 

5.14 

1507 

4 

Core  Volume  3 

7.3 

5.88 

1507 

5 

Core  Volume  4 

7.3 

6.62 

1506 

6 

Reactor  Outlet  Plenum 

74 

7.36 

1502 

7 

Hot  Leg  Piping 

19 

7.56 

1500 

8 

Pressurizer 

138 

6.24 

1500 

9 

Surge  Line 

4.9 

7.78 

1500 

10 

Steam  Gen.  Inlet 

18 

7.63 

1498 

11 

Tube  Volume  1 

5.1 

10.21 

1497 

12 

Tube  Volume  2 

5.1 

12.34 

1496 

13 

Tube  Volume  3 

5.1 

14.47 

1495 

14 

Tube  bends  4 

7.3 

16.60 

1493 

15 

Tube  Volume  5 

5.1 

14.47 

1490 

16 

Tube  Volume  6 

5.1 

12.34 

1489 

17 

Tube  Volume  7 

5.1 

10.21 

1488 

18 

Steam  Gen.  Outlet 

15 

7.63 

1485 

19 

Cold  Leg  Piping 

16 

7.63 

1516 

20 

Reactor  Downcomer 

44 

0.0 

1508 

21 

Primary  Containment 

23016 

0.0 

15 

*Dimension  from  reference  level  to  lower  elevation  of  volume 
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TABLE  B-l  (coat’d) 


SUMMARY  OF  RELAP-3  INPUT  DATA  FOR  MH-1A  MODEL 
SINGLE  ENDED  HOT  LEG  BREAK 


Junction  Data 


Minimum  Junction 

Junction  From  To  Flow  Area  Height 

No.  Volume  Volume _ (ft2) _ (ft) 


1 

1 

2 

9.67 

4.4 

2 

2 

3 

9.67 

5.2 

3 

3 

4 

9.67 

5.9 

4 

4 

5 

9.67 

6.6 

5 

5 

6 

9.67 

7.4 

6 

6 

7 

0.667 

9.1 

7 

8 

9 

0.18 

8.2 

8 

9 

7 

0.18 

8.2 

9 

7 

10 

0.667 

8.0 

10 

10 

11 

2.38 

10.2 

11 

11 

12 

2.38 

12.3 

12 

12 

13 

2.38 

14.5 

13 

13 

14 

2.38 

16.6 

14 

14 

15 

2.38 

16.6 

15 

15 

16 

2.38 

14.5 

16 

16 

17 

2.38 

12.3 

17 

17 

18 

2.38 

10.2 

18 

18 

19 

.473 

8.1 

19 

19 

20 

.667 

9.1 

20 

20 

1 

1.046 

0.32 

21 

(Fill  System) 

20 

0.0089 

11.6 

22 

7 

21 

0.667 

9.5 

Initial  Parameters 

Power 

49.5 

MWt 

Coolant  heating  rate 

1.69xl08 

Btu/hr 

Flow  rate 

1080 

lb  /sec 

Total  mass 

18,407 

lb 

Total  energy 

1.03xl07 

Btu 
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APPENDIX  C 


PUMP  REPRESENTATION 


The  effects  of  the  reactcr  coolant  pimps  on  the  system  blowdown  are 
represented  In  the  analysis.  The  primary  effects  simulated  include  the 
pump  head  versus  flow  characteristics,  punp  cavitation,  coastdown  fol¬ 
lowing  pump  trip,  and  flow  resistance  across  the  pumps  after  shutdown. 
The  reactor  coolant  pusps  in  the  MH-1A  consist  of  two  single  stage 
centrifugal  punps  connected  in  parallel  by  a  piping  volute  assent) ly. 
Table  C-l  shows  the  nominal  speed  pimp  head  versus  flow  curve  for  paral¬ 
lel  performance  used  to  represent  the  primary  pumps  in  RELAP-3.  This 
curve  is  based  on  the  manufacturer's  test  data  (reference  3)  conserva¬ 
tively  extended  to  the  range  of  flows  encountered  in  a  LOCA. 


The  loss  of  effectiveness  of  the  punps  due  to  cavitation  and  flow 
coastdown  (if  pump  trip  is  assumed)  is  evaluated  in  RELAP-3  from  the 
following  equation: 


Ap  Pu“P  =  144 


Ap  pump  ^l  •  F2 


(C-l) 


where: 


pump 

PA 

fl 

F1 

PD 


PV 

F2 

*2 


normal  punp  head  from  flow  characteristic  curve,  ft 

fluid  density  at  pump  junction,  lb /ft3 

1.0  -  Fi  (PD-Pi)2,  loss  of  effectiveness  0<fi <1.0 

cavitation  constant  (enpirical) 

pressure  at  inception  of  cavitation 

Pv  +  F2p* /144 

vapor  pressure  at  pump  inlet,  psia 
NPSH,  ft 

pump  coastdown  multiplier 


The  term  f2  simulates  the  loss  of  effectiveness  due  to  pump  coastdown  and 
is  applied  if  loss  of  power  to  the  pumps  is  assumed.  If  no  pump  trip 
occurs  then  f2  remains  equal  to  one.  The  coastdown  multiplier  for  loss 
of  power  to  both  pumps  was  developed  from  an  approximation  of  the  pump 
characteristics  and  is  consistent  with  flew  coastdown  reported  in 
Reference  11.  The  shutdown  resistance  of  the  pumps  following  flow 
coastdown  was  determined  from  manufacturer's  test  data  and  is  defined 
in  RELAP-3  as: 
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where: 

Ap  pump  =  pump  pressure  differential,  psi 
PL  **  coolant  density  at  punf>  junction,  lb/ft3 

tfi  =  average  flow  across  pump  junction,  lb/sec 

kg  *  punj>  shutdown  resistance 

Table  C-2  shows  calculated  shutdown  resistances  based  on  the  above  equation 
for  various  pump  impeller  conditions.  Also  shown  are  additional  parameters 
used  to  represent  the  MH-1A  coolant  pumps  in  the  analysis.  The  pump  net 
positive  suction  head  (NPSH)  specified  by  the  manufacturer  was  assumed  and 
a  conservative  cavitation  constant  was  selected  so  as  to  result  in  a  rapid 
loss  of  effectiveness  upon  cavitation. 

TABLE  C-l 

MH-1A  NOMINAL  SPEED  PUMP  HEAD  VERSUS  FLOW  INPUT  DATA 
(Parallel  Performance) 


Flow 

gpm 

Head 

ft 

-8000 

170 

+2000 

160 

4000 

152 

6000 

144 

8000 

132 

10,200 

110 

12,000 

80 

14,000 

44 

, 

16,000 

0 

100,000 

0 
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HH-iA  PRIMARY'  COOLANT  PUMP  DATA 


Number  of  Pumps 


2 


Type 


Centrifugal 


Pump  NPSH,  ft 


100 


Cavitation  constant  (empirical) 


.0018 


Pump  shutdown  resistance*  (490°F) 
Free  rotors 
One  locked  rotor 
Two  locked  rotors 


5.0  x  10-4 
7.9  x  10"4 
10.8  x  10“* 


*From  equation  C-2  based  on  manufacturer's  test  data. 
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REACTOR  KINETICS 


The  reactor  kinetics  calculation  in  RELAP-3  is  a  point  kinetics  model 
with  appropriate  reactivity  feedback  effects.  The  standard  reactor  kinetics 
equations  used  in  the  program  are: 


6 


dn 

e*  [£.  -  i]  n 

dt  " 

i*  $* 

* 

dCi  „ 

-Ai  C±  +  Bi 

.  n 

dt 

i* 

r 

i=l 


Ai  C± 


where: 

t  »  time,  sec 

n  *  reactor  fission  power,  Btu/sec 

p  »  total  reactivity,  Ak/k 

*  effective  delayed  neutron  fraction  of  the  i*1*1  group 
0*  a  effective  delayed  neutron  fraction 
*i  *  decay  constant  of  i**1  delayed  neutron  group,  sec"1 
l*  ■  neutron  generation  time,  sec 

Gj_  «  precursors  concentration  for  i^1  delayed  neutron  group 

Also  included  is  a  determination  of  gamma  heating  due  to  decay  of  the 
fission  products.  The  total  power  is  calculated  as  the  sum  of  the  contri¬ 
butions  from  the  instantaneous  fission  rate  and  decay  of  the  fission 
products.  All  power  is  assumed  to  be  generated  in  the  fuel  elements, 
direct  gamma  heating  of  the  coolant  not  being  considered. 


fhe  major  reactivity  feedback  effects  included  In  the  kinetics  equa¬ 
tions  are  contributions  from  control  rod  motion,  moderator  voids,  moderator 
temperature,  and  Doppler  broadening.  The  reactivity  is  input  either 
explicitly  as  a  function  of  time  as  with  reactivity  insertion  due  to  reactor 
trip  or  implicitly  through  core  feedback  mechanisms  determined  from  coolanc 
density,  coolant  temperature,  and  fuel  temperature  for  each  core  volume. 

The  total  reactivity  is  determined  by  summing  up  the  individual  contributions. 
The  initial  value  of  reactivity  is  assumed  to  be  zero  corresponding  to  a 
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steady-state  equilibrium  condition.  Reactivity  feedback  from  the  moderator 
teaq>erature  and  voids  is  expressed  in  terms  of  a  linear  negative  feedback 
coefficient  clossn  so  as  to  result  in  a  conservative  power  transient  during 
blowdown.  Nonlinear  Doppler  feedback  is  input  in  the  form  of  a  curve  of 
reactivity  versus  effective  fuel  element  temperature  based  on  the  Doppler 
coefficient  shown  in  Figure  D-l.  Table  D-l  summarizes  the  reactivity  coef¬ 
ficients  and  nuclear  constants  which  were  used  to  represent  the  MH-1A  in 
the  reactor  kinetics  calculations  in  RELAP-3.  The  dominant  influence  on 
the  core  pcwer  transient  results  from  the  moderator  void  feedback  coeffi- 
cient  and  this  coefficient  is  input  as  a  conservatively  low  negative  value. 


TABLE  D-l 


REACTOR  KINETICS  DATA  FOR  MH-1A  CORE 


£  Mean  neutron  lifetime,  sec 
6*  Effective  delayed  neutron  fraction 

61  'l 

62  ; 

63  k  Effective  Neutron  Fraction  for  i“*  Group 

64 

65  1 

66  j 


A3  ‘>  Decay  Constant  for  i^1  Group,  sec 
A  4  : 

A5 


aT  Moderator  Coefficient  (BOL)  dk/°F 

zero  power  (68°F) 
full  power  (490° F) 

ajj  Doppler  Coefficient  (See  Fig.  D-l)  dk/°F  (UO2) 

zero  power  (68°F,  0  MW) 
mil  power  (4903F,  45  MW) 

ay  (1)  Void  Coefficient  dk/Z  void 


1.3  x  10~5 
.007015 

.00027 

.00149 

.00132 

.00285 

.00090 

.00018 

.0127 

.0317 

.115 

.311 

1.40 

3.87 


-0.6  x  10~4 
-2.5  x  10*4 


-2.6  x  10'5 
-1.6  x  10~5 

-2.2  x  10'3 


(1)  Nominal  value  -  reduced  by  25%  in  analysis 
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FIGURE  D-l  DOPPLER  COEFFICIENT  VS  EFFECTIVE  FUEL  TEMP K RAT 


APPENDIX  E 


TIME  STEP  SIZE 


RELA°-3  utilizes  a  forward  finite  difference  approximation  as  its 
numerical  advancement  technique.  Since  the  same  time  step  size  is  used 
for  all  volumes  and  junctions,  the  stability  of  the  solution  is  controlled 
by  the  most  sensitive  element  of  the  model.  The  time  step  size  must  be 
smaller  than  the  natural  oscillation  period  of  this  element  by  a  factor 
of  5  to  10  or  numeric  instabilities  can  result.  An  appropriate  time  step 
is  given  by: 

A  t  *(2ir/na)i/  V  £/A 
where 

n  «  5  to  10 

a  *=  Velocity  of  sound  in  medium  contained  in  volume 

V  “  Volume  of  fluid 
i/A  *  Path  inertia 

In  the  analysis  the  most  restrictive  element  was  determined  for  each 
model  and  the  time  step  for  stability  was  calculated  using  the  element 
dimensions  and  the  acoustic  velocity  of  the  fluid  in  the  element.  A 
factor  of  n  **  10  was  used.  The  accuracy  of  the  computer  solution  was 
verified  by  reducing  the  time  step  by  one  half  and  observing  the  calculated 
results.  Comparison  of  the  initial  core  flow  rates  and  fuel  temperatures 
showed  no  significant  differences.  Increasing  the  time  step  above  the 
criterion  resulted  in  numerical  instability.  Other  cases  used  initial 
time  steps  consistent  with  the  stability  criterion  with  the  factor  n  -  10 
to  reduce  the  time  step  below  the  natural  frequency. 

Later  in  the  problem  solution,  ether  types  of  numeric  instability 
due  to  mass  and  energy  flow  could  occur.  The  time  steps  in  this  region 
of  the  problem  were  selected  so  as  to  prevent  the  fluid  mass  or  energy 
transfer  during  a  time  increment  from  becoming  large  relative  to  mass 
or  energy  contents  of  the  element. 
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APPENDIX  F 


FUEL  ROD  PROPERTIES 

The  variation  of  fuel  rod  properties  with  temperature  was  considered 
in  the  analysis.  The  high  temperature  property  data  incorporated  in 
THETA  1-B  is  summarized  in  Table  F-l  and  Figures  F-l  and  F-2.  The  figures 
show  the  thermal  conductivity  and  heat  capacity  of  stainless  steel  and 
UO2  fuel,  respectively.  Data  for  stainless  steel  was  taken  from 
Reference  23.  The  temperature  dependent  conductivity  was  based  on 
composite  data  for  austenitic  stainless  steel  with  the  error  considered 
+  5  percent.  Above  1330°F  the  thermal  conductivity  was  conservatively 
assumed  constant  since  the  rate  of  increase  may  be  limited  by  oxide 
formation.  The  heat  capacity  of  the  cladding  was  computed  based  on  the 
chemical  composition  of  the  alloy  assuming  18  percent  chromium  and  8 
percent  nickel  with  the  remainder  being  iron.  Heats  of  transformation 
corresponding  to  phase  changes  in  Fe  were  assumed  at  temperatures  of 
approximately  1400  and  1670°F.  The  melting  temperature  of  the  cladding 
was  taken  as  2550°F,  the  recommended  lower  limit  for  austenitic  stainless 
steel.  The  UO2  property  data  in  THETA  1-B  was  based  on  the  information 
in  Reference  23.  The  correlations  used  in  the  code  are  reported  in 
Reference  8. 


TABLE  F-l 


HIGH  TEMPERATURE  PROPERTY  DATA  OF  FUEL  ELEMENT  MATERIALS 


Stainless  Steel 
(Ref.  23) 

U02 

(Ref.  8.23) 

1.  Melting  point,  °F 

2552-2659 

5072 

2.  Latent  heat .of  fusion 
Btu/lb 

Cr  Fe  Ni 

121  119  129 

67.5  +  1.8 

3.  Coefficient  of  thermal 
expansion  in/in  °F 

1.12  x  10"5 
mean  coefficient 
(32  -  1832°F) 

3.7?.  x  10"5+1.787  x 
10-9  T  Tin  °F 
(32-4082°F) 

4.  Density  @  77°F, 
lb /ft 3 

493 

(95% 

650 

theoretical  density) 

5.  Thermal  conductivity 
Btu/hr  f t°F 

<32-1364°F) 

Fig-  F-l 

Fig.  F-2 

6.  Heat  capacity,  Btu/lb°F 

Fig.  F-l 

Fig.  F-2 
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Heat  Capacity,  Btu/lb  °F  Thermal  Conductivity,  Btu/hr-ft  °F 
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Temperature  »°F 


Stainless  Steel 
Heat  Capacity 


500  1000  1500  2000  2500 

Temperature,  °F 

FIGURE  F-l  STAINI.ESS  STEEL  CC  CTIVITY  AND  HEAT  CAPACITY 

VS  TEMPER/ ; JAE 
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FIGURE  F-2  U02  CONDUCTIVITY  AND  HEAT  CAPACITY 
VS  TEMPERATURE 
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appendix  g 


METAL  WATER  REACTION 


When  raised  to  high  temperatures  in  the  presence  of  steam  or  water 
vapor,  stainless  steel  cladding  can  react  with  steam  to  form  metal  oxides 
and  hydrogen.  This  reaction  is  exothermic  and  the  energy  release  can 
contribute  to  the  clad  temperature  transient.  The  effects  of  the  stain¬ 
less  steel  -  water  reaction  are  represented  in  THETA1-B  by  means  of  the 
parabolic  rate  law.  The  reaction  between  stainless  steel  and  steam 
proceeds  relatively  slowly  at  low  temperatures,  but  the  rate  rapidly 
increases  as  metal  approaches  its  melting  temperature.  The  recommended 
equation  for  stainless  steel  (Reference  23)  In  terms  of  the  weight  gain 
due  to  oxidation  of  the  clad  above  180CaF  is: 

w2/t  **  2.4  x  1012  exp  [84,300  2400] ,  mg2  (0)/cm4sec  (G-l) 

RT 


where  w  is  the  weight  gain  (oxygen)  per  unit  surface  area  in  mg/ cm2, 
t  is  the  time  in  seconds,  R  is  the  gas  constant  in  cal /mole  °K,  and 
T  is  tue  absolute  temperature.  The  heat  of  reaction  for  the  above 
process  is  about  254  cal/gm  of  stainless  steel  oxidized  in  the  solid 
state  (Reference  24).  For  purposes  of  analysis  the  major  reaction  pro¬ 
duct  formed  is  assumed  to  be  Fe^1*  although  other  oxides  may  be  present. 
In  THETA1-B  the  rate  equation  is  expressed  in  terms  of  the  amount  metal 
reacted  as  follows: 


dr 

dt 


_ 5 -  exp  (-AE/RT) ,  in/sec 

<r-rQ) 


(G-2) 


where: 


r  *  radius  of  unreacted  metal  in  clad,  in 
r0  -  original  radius  of  un reacted  clad,  in 
t  »  time,  seconds 
K  *  rate  law  constant,  in2 /sec 
AE  *  activation  energy,  cal/mole 
R  “  gas  constant,  cal/mole  °K 
T  *  cladding  temperature,  °K 


In  the  above  relation,  (r-rQ)  represents  the  oxide 
layer  which  retards  the  reaction  in  proportion  to  its  thickness.  The 
rate  constant  (K)  is  related  to  the  pre-expotential  factor  of  the  rate 
equation  (A)  as  follows : 


2pm*  (2.54)* 


in2 /sec 


(G-3) 
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where  A  is  expressed  in  terms  cf  mass  of  metal  (Fe)  reacted  per  unit  surface 
area  per  sec  quantit;  squared  in  gn2/cm4sec2  and  pm  is  the  density  of  stain¬ 
less  steel  in  gm/cc.  At  relatively  high  temperatures  the  reaction  may  be 
limited  by  the  availability  of  steam  at  the  oxide  sutface.  For  conserva¬ 
tism  no  steam  limiting  is  permitted,  in  the  analysis,  i.e.  ,  the  rate  of 
reaction  is  assumed  to  equal  the  rate  calculated  by  equation  G-2.  At 
temperatures  approaching  the  melting  point  of  stainless  steel  0v2550°F) 
a  sudden  rapid  increase  in  surface  area  and  reaction  rate  is  observed. 

For  engineering  purposes  it  is  conservatively  assumed  that  100  percent 
reaction  occurs  if  the  cladding  reaches  the  melting  temperature.  Actually 
the  percent  reaction  would  again  be  limited  by  diffusion  of  steam  to  the 
reacting  surface.  Table  G-i  summarizes  the  metal  water  reaction  constants 
used  in  the  analysis. 


TABLE  G-l 

METAL  WATER  REACTION  CONSTANTS 


Material 

Melting  point,  °F 

Density,  gm/cc 

Heat  of  reaction,  cal/gm 

Ratio  of  mass  of  Fe  to 

mass  of  oxygen  reacted 

AE  activation  energy,  cal/mole 

R  rate  law  constant,  in2/sec 

R  gas  constant,  cal/mole  °K 


Stainless  Steel 

2550 

8 

254 

2.61 

84,300  +  2400 
1.985  x  104 
1.987 
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